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Abstract 

A growing part of the behavioral finance literature has addressed some of the stylized facts of 
financial time series as macroscopic patterns emerging from herding interactions among groups 
of agents with heterogeneous trading strategies and a limited rationality. We extend a stochastic 
herding formalism introduced for the modeling of decision making among financial agents, in 
order to take also into account an external influence. In particular, we study the amplification 
of an external signal imposed upon the agents by a mechanism of resonance. This signal can be 
interpreted as an advertising or a public perception in favor or against one of the two possible 
trading behaviors, thus periodically breaking the symmetry of the system and acting as a contin- 
uously varying exogenous shock. The conditions for the ensemble of agents to more accurately 
follow the periodicity of the signal are studied, finding a maximum in the response of the system 
for a given range of values of both the noise and the frequency of the input signal. 
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1. Introduction 



The ubiquitous observation of stylized facts or non-trivial universal statistical properties 



(IConn. 12001) in financial time series has traditionally been explained within the efficient mar- 
ket hypothesis as a direct consequence of the statistical properties of the news arrival process. 
Recent years, however, have witnessed the appe arance of a growing number of contributions 
based on heterogeneous interacting agents (.Hommes , 20061) interpreting these stylized facts as 
the macroscopic outcome of the diversity among the economic actors, and the interplay and 
connections between them ( Kirman Ul992 : Lux and Westerhoff , 20091: Farmer and Geanakoplos 
2009t IColander et al.L i2009). Heterogeneity refers here to the agents' capability to choose from a 
set of different market strategies or trading rules. Regarding their interplay, different interaction 
mechanisms have been studied, whether direct or indirect, global or local. 

We will focus on a series of stochastic models of information transmission inspired by 
KirmanI (ll99lUl993h . and whose main ingredient is their emphasis on the processes of social 
interaction among agents, based on herding behavior or a tendency to follow the crowd. These 



* Corresponding author Tel.: +34 971 25 98 82; Fax: +34 971 17 32 48 
Email address: adrian.carroiaifisc.uib-csic.es (Adrian Carro) 



Preprint submitted to arXiv 



February 27, 2013 



mechanisms of social interplay have also dra wn a growing attention from different disciplines 



during recent decades (ICastellano et al.l 120091) 



Inspired by a series of entomological experiments with ant colonies, Kirr nan proposed a 
stochastic herding formalism to model decision making among financial agents ( lKirmanlll993h . 
In the experiments with ants, entomologists observed some asymmetric macroscopic patterns 
emerging from an apparently symmetric situation: when ants were faced with a choice between 
two identical food sources, a majority of the population tended to exploit only one of them at a 
given time, switching its foraging attention focus to the other source every once in a while. In 
order to explain this behavior, Kirman developed a stochastic model where the probability for an 
ant to change its foraging source is a combination of two terms, one related to pairwise herding 
interactions or recruitment probability, and the other as an autonomous switching tendency or 
idiosyncratic behavior This very simple model can also be interpreted in terms of market be- 
havior by just replacing an ant's binary choice between food sources by a market agent's choice 
between two trading strategies. 



A series of subsequent papers (ILux and Marchesil 119991: lAlfarano et al.L I2005L 120081 1201 Ih 



has focused on explaining some of the stylized facts observed in empirical data from finan- 
cial markets in terms of herding models of the Kirman type. Some of the universal features 
of financial time series which have been reproduced in this framework are: the non-Gaussian 
or leptokurtic character of the unconditional distribution of returns — also known as "fat tail" 
shape — , the temporal bursting behavior of the volatility — or volatility clustering — , and the 
positive autocorrelation of absolute and squared returns. 

The main purpose of the present paper is to contribute to the theoretical assessment of the 
influence of an external significant agent, media group or institution upon a financial market. To 
this end, we study an agent-based herding model of the Kirman type extended to take into account 
incoming information. In particular, we add an external periodic signal which will, at regular 
time intervals, induce a modification of the relative attractiveness of both possible options, thus 
helping the agents to switch their choice in one direction or the other alternatively. This external 
forcing can be interpreted, in terms of financial markets, as an advertising or a public perception 
in favor or against one of the two possible trading behaviors, thus periodically breaking the 
symmetry of the system. The conditions for the ensemble of agents to more accurately follow 
the periodicity of the signal are studied. In particular, we measure the response of the system as a 
function of the tendency of the agents to autonomously switch their trading strategy. We find, for 
a wide range of the model parameters, a resonance phenomenon, i.e., a maximal amplification 
of the incoming information as a function of a given parameter of the model. Even if the perfect 
periodicity of the external signal may not be a realistic feature of real financial markets, it does 
not affect the generality of our results concerning the conditions for the system to better foUow 
any given signal. 

The outline of the paper is as follows. In section |2] we present the original Kirman model, 
derive the Langevin equation describing the process, find its stationary solution and describe 
the dynamics in terms of an effective potential. The motivation for a Langevin description is to 
clearly observe the role played by the parameters of the model, as well as for a direct comparison 
with the externally forced case, which is analytically developed in section [3] The numerical 
results obtained with this latter model, their analysis as well as their graphical presentation, are 
covered in section|4] where we also discuss some particular features of the resonance found. We 
draw, in section |5] the main conclusions and add some general remarks. 

It is also worth mentioning that the method we used to solve the model is different from those 
used in the previous hterature. We apphed a .Gillespie, (.1977,) algorithm in order to generate 
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statistically correct trajectories of the stochastic equation. Therefore, we find the unbiased times 
when the transitions take place. A brief description of this simulation method will be left for a 
final appendix. 



2. The original herding model 

A stochastic model of information transmission was proposed rrmani (Il993h in order to 
explain the collective behavior observed in ant colonies when they faced a choice between two 
identical food sources located in the neighborhood of their nest. Despite an apparently symmet- 
ric situation, being both sources exactly equal and located at the same distance, entomologists 
observed the systematic emergence of asymmetric macroscopic patterns of foraging among the 
ants. Indeed, at any point in time, a clear majority of the ant population is found exploiting one 
of the food sources, while a minority concentrates on the other on^ Moreover, the experiments 
showed that every once in a while a switch occurs in the system and a majority is found foraging 
in the source previously neglected. Therefore, averaging over time, the probability distribution 
of the fraction of ants visiting one or the other food sources becomes bimodal. 

In order to explain this behavior, Kirman postulates, on the one hand, the existence of a herd- 
ing propensity among the ants, i.e., a tendency to follow the crowd, which implies the existence 
of some kind of interaction among them with information transmission. This communication 
process has been found by entomologists and works mainly through hormone signals or physical 
recruitment. On the other hand, he also assumes the ants to randomly explore their neighbor- 
hood looking for new food sources, so every one of them has an autonomous or idiosyncratic 
switching probability due to this stochastic search. Thus, the model developed by Kirman basi- 
cally states that the probability for an ant to change its foraging source is a combination of two 
terms, one related to pairwise herding interactions or recruitment probability, and the other as an 
idiosyncratic switching tendency. 

This model can be easily interpreted in terms of market behavior and has been used to 
mod el contagion and herding phenomena in financial markets in a number of subsequent pa- 
pers (lAlfarano et all l2005l l2008l lAlfarano and MilakovicI EoOOl lAlfarano et all l201lb . some of 
them including significant variations (|Lux and Marchesi, 1999). In a market context, an ant's 
binary choice between food sources is replaced by a market agent's choice between two trading 
strategic^ For instance, in a foreign exchange market, agents can adopt different trading strate- 
gies such as a fundamentalist or a chartist forecast of future exchange rate movements. A further 
example would be the choice between an optimistic or pessimistic tendency among the chartist 
traders . In these examples, the Kirman model would be the decision making mechanism among 
financial agents, who decide whether to buy or sell in a given situation, thus giving rise to market 
switches between a dominance of one or the other strategy. 

Let us now briefly review the formalization of the original Kir man model and the de riva- 
tion of some analytical results a long the lines of previous works by lAlfarano et al.l (l2008b and 
Kononovicius and GontisI (12012 V Let the market, or ant colony, be populated by a fixed number 
of agents A^, and let n be the number of those agents choosing one option, while N — n choose the 
other one. For the sake of clarity, we will refer hereafter to the agents n as optimistic and to the 



' So as to make sure that the initial situation is symmetric, the entomologists used also a setting with only one food 
source and two identical and separated paths, obtaining equivalent results. 

^These trading strategies may be related to some particular rules for the formation of his expectations about the future 
evolution of prices, or result from differences in access and interpretation of present and past information. 
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agents N — n as pessimistic. Thereby, n £ 0, 1 , ...,A^ defines the configuration of the system. The 
evolution of the system is given by two terms: on the one hand there are pairwise encounters of 
agentfl after which one of them may copy the strategy of the other, and on the other hand there 
are idiosyncratic random changes of state. 

There have been mainly two different implementations of the pairwise encounter term in 
the literature. In his seminal 1993 paper, Kirman proposed the use of a herding intensity that, 
for each agent, is proportional to the fraction of agents in the opposite state. One of the main 
drawbacks of this original formalization has been pointed out to be its lack of robustness with 
respect to an enlargement of the system size, or N-dependence, since an increasing number of 
participants in the market causes the stochasti city to vanish and th erefore th e stylized facts to fade 
awa> 1. On the contrary, som e later authors ( Alfarano et al.i [2005 1 12008: Alfarano and Milakovic , 



2009tlAlfarano et al.Ll201 11) avoided this problem with an alternative modelization of the pairwise 
interaction mechanism with a herding intensity that, for each agent, is proportional to the total 
number of agents in the opposite state, thus allowing each individual to interact with any other 
regardless of the system size. We will hereafter adopt the second and more recent formalism. 

An additional assumption of the model is the lack of memory of the agents, so their prob- 
ability of changing state does not depend on the outcome of previous encounters, neither on 
previous idiosyncratic switches. Therefore, the stochastic evolution of the system can be for- 
malized as a Markov process depending, at each time step, just on the probability to switch 
from the present configuration of the system n to some other state n' in a time interval At, 
denoted by P{n' ,t + At\n,t). However, if this time interval At is taken to be small enough, 
then the probability to observe multiple jumps is negligible and we can restrict our analysis 
to n' = n ± 1 . Furthermore, the probabilities would then be related to the transition rates per 
unit time as P{n' ,t + At\n,t) — n{n — > n')At. According to the evolution of the model described 
above, the transition rates for each individual /, n'^ — Kiipessimistic optimistic) and — 
ni{optimistic pessimistic), can be formally defined as 

= a + hn, 

%^ =a + h{N-n), 

where the parameters a and h stand for the idiosyncratic switch and the herding intensity coeffi- 
cients respectiveljQ. The rates for the whole system are, therefore, 

7r+(n) = 7r(n— s>n + l)= {N-n){a + hn), 

(1) 

n («) = n{n^n-\)= n{a + h{N-n)). 

There are two parameters in the model, a and h, but one of them can be used as a rescaling of the 
time variable, so that there is only one relevant parameter, such as £ = |. 



^Having one-to-one meetings implicitly assumes the density of agents to be low enough so as to avoid multi-agent 
meetings. 

^By implem enting the dynamics in different network structures and hierarchies, Alfarano and Milakovi^ J2009ft ; 
lAlfarano et alj 420 1 1 ) pointed out that the problem of N-dependence is related to the structural heterogeneity of the 
model. 

^Obviously, only one of these transition rates would apply at each time step for each agent in case agent i is 

currently in the pessimistic or — 1 state, and Kj' if it is currently in the optimistic or + 1 state. 
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Using the one-step increment and decrement operators!!, E and E"' (TVan KampenL l2007l) . 
and calling P{n,t) the probability to have n agents in the optimistic state at time t, the master 
equation describing the process is written as 



dP{n,t) 
dt 



= (E - 1) [K-{n)P{n,t)\ + (E-i - 1) [K+{n)P{n,t)\ . 



(2) 



We now change from the extensive variable n in the range n E [0,N] to an intensive variable 
m — 2fj — l, giving the population configuration or opinion index in the range m G [— 1 , + 1 ] . This 
intensive variable m can be thought of as continuous if the number of agents is large enough. 
Note, however, that the limit of an infinite number of agents is never the case in real soc ial and 
economical systems, where finite size effects may play a role (^ Toral and Tessone , 2007 ). Note 
as well that a population configuration m = would imply a perfect balance of opinions, while 
m — — l and m = +l would signal a total agreement on the pessimistic and the optimistic opinions 
respectively. Forthis new variable, the partial derivation is ^ = (^) ^ and the relation between 
the old and the new variable probabilities becomes P{n,t) — P{m,t)^. If we also take into 

account that the one-step increment operator can be expressed as E = exp | ^ | and, furthermore, 
we assume a large number of market participants A^, then this one-step increment operator term 
can be approximated by a Taylor series up to second order in ^ as 
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Introducing this and the transition rates dl} into the expression (|2]), the master equation can be 
approximated by a Fokker-Planck equation. 



dP{m,t) 
dt 



d_ 

dm 



2amP{m,t] 



1 d^ 



2dn 



4a 

— +2h{\-m'-) ]P{m,t) 



(3) 



where /x(m) = —2am plays the role of a drift term and D{m) = ^+2h{l —m^) is the diffusion 
coefficient. 

The steady state solution of the Fokker-Planck equation (O is 



P,,{m) - 3f- 



a 
2Nh 



(l-m2) 



(4) 



where iF^' is a normaUzation factor. Observing the functional form of this solution, one can 
notice that the sign of the exponent will determine whether the probability distribution is uni- 
modal with a peak centered at m = or bimodal with peaks at the extremal values m — —l and 
m = + 1 . Therefore, when the idiosyncratic switching a is larger than the herding intensity h, we 
find a unimodal distribution, meaning that, at any point in time, the most likely outcome of an 
observation is to find the community of agents equally split between both options. On the con- 
trary, when the herding h exceeds the idiosyncratic switching intensity a, a bimodal distribution 
is found, meaning that, at any point in time, the most likely outcome of a static observation is to 



^The increment operator E is defined as £[/(«)] = f{n +1), wliere f{n) is an arbitrary function of the integer argu- 
ment n, whereas the decrement operator E"' has the effect E^' [/(h)] = f{n — 1). 
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find a large majority of agents choosing the same option. Nevertheless, in different observations, 
the option chosen by the majority may be different. Note as well that when a — h the probability 
distribution is uniform. Because of the ergodicity of the model, these probability distributions 
can also be understood in terms of the fractional time spent by the system in each state. 

For a large number of agents A^, we can neglect the first term inside the brackets in the steady 
state solution (|4|l, and thereby write the normalization factor a^ 

1 ^ r{2a/h) ^ r(2e) 

r(fl//!)2 r(e)2 ■ 

Returning to the Fokker-Planck equat ion Q and applying the usual transformation rule 
( Van Kampenl 2007 ). within the It^ 1 1951 ) convention, we find the Langevin or stochastic dif- 



ferential equation describing the process. 



til = jU (m) + y/D{mj • ^ (f ) 

f4-a 

^-2am+ a/ — + 2/i(l-w2)-^(f), 

where £, (?) is a Gaussian white noise, i.e., a random variable with a zero mean Gaussian distri- 
bution, (f )) = 0, and no autocorrelation, (^ (r)^ (f')) = 5(f - ?')• 

Let us first analyze the role played by the noise or diffusive function inside the square root, 
D{m) — ^ +2/i(l — m^). The first term, dependent on a and inversely proportional to the sys- 
tem size A^, is related to the "granularity" of the system, and so it vanishes in the continuous 
limit N ^ °°. It basically states that for any finite system there are always finite-size stochastic 
fluctuations related to the fact that the agents have the ability to randomly change their choice. 
The second term of the diffusion function is a multiplicative noise term, i.e., a noise whose in- 
tensity depends on the state variable itself. Furthermore, it is the only term dependent on the 
herding coefficient, so we will refer to it hereafter as herding term. As this multiplicative noise is 
maximum for m — and vanishes for ;« = ±1, it tends to move the system away from the center 
and towards those extremes by making a random partial agreement on one or the other possible 
opinions grow to a complete consensus. Note that for a — the consensual states m ~ ±1 be- 
come absorbing states. We will nonetheless refer to them as absorbing states also for a 0, even 
if they are not properly absorbing in that case. 

Regarding the deterministic drift term, ji (m) ~ ~2am, it drives the system back to a balanced 
position at the center of the opinion index, m = 0. Therefore, we have a competition between 
two driving forces. One of them is of a stochastic nature, it is dominated by the herding term 



'Certainly we can also apply the probability nomialization condition without neglecting the N-dependent term. In 
this last case, the normalization factor iF^ ' would be 
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where B[x\a,b\ is the incomplete beta function, given by 
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for large and it tends to favor the formation of a majority of agents sharing the same opinion. 
Whereas the other is of a deterministic nature, it is related to the idiosyncratic switches and 
it tends to break these majorities and drive the system back to a balanced situation, where the 
agents are equally distributed between both opinions. Depending on the relative magnitude of the 
idiosyncratic and the herding parameters, a and h, one or the other behavior prevails. Indeed, the 
particular functional form of the noise in the Kirman model induces a transition in the dynamics 
of the system, from a monostable to a bistable behavior when increasing the value of h relative 
to that of a. 

In order to clearly observe this noise-induced transition, let us directly write here an effective 
potential for the Fokker-Planck equation Q which takes into account the effect of the multiplica- 
tive noise over the deterministic potential related to the drift term, 

V^fi(m) = {h - a) ln(l - m^). (6) 

The derivation of this effective potential, basically consisting of assuming an exponential func- 
tional form for the stationary probability distribution, is left for [Appendix B| It is obvious, from 
equation Q and Figure [T] that a transition will occur at h — a from a one well to a double well 
potential, for increasing h. 

Effective Potential Vg^m) 
a < h a = h a > h 



its 



0.00 




1-1 1-1 1 

m m m 

Figure 1: Effective potential Veff(m) for a herding intensity ho = 0.01 and three different values of the idiosyncratic 
switching tendency, a = 10"' {a <h), a = 10"^ (a = h) and a = 10"' {a > h). 

Another interesting feature worthy of remark here is the fact that, although both wells at the 
extremes of the opinion space in the effective potential are theoretically infinite, in fact, they are 
not absorbing states and so the system will leave them with rate proportional to a. This can be 
understood by a reexamination of the Langevin equation (|5), where we notice that precisely at 
the extremal states m ± 1 the only term acting upon the system is the deterministic force towards 
the center of the opinion space. To sum up, the role of the herding h is to induce a bistable 
effective potential, while that of the idiosyncratic switching a is to allow for transitions between 
both extremal wells. 

An example of a stochastic realization of equation (|5]l is shown in the panel a of Figure [3] 
for a herding intensity ten times larger than the idiosyncratic switching coefficient, £ = f = jq- 
We can observe, in this panel, the tendency of the system to be temporarily absorbed around the 
consensual states with random switches between them. Some of those switches are not successful 
and the system returns to the previous consensus before reaching the opposite one. 
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3. The model with an external signal 



We now proceed to extend the model in order to account for an input of external information 
able to influence the choices made by the economic actors. We are concerned here with the 
search for the conditions which allow the market to more accurately follow an information arrival 
process. To this end, we basically add a weak periodic input signal or driving force, which will, 
at regular time intervals, induce a modification of the transition rates helping the switches in 
one or the other direction alternatively. The only reason for the signal to be perfectly periodic 
is to simplify the measurement and quantification of the similarity between the input driving 
signal and the output results of the model in terms of market behavior. Nevertheless, our results 
concerning the search for the parameter ranges for the market to more accurately follow the 
external signal are general and independent of its particular shape. The only relevant features of 
the driving signal which have a significant effect on the results are its strength and its frequency 
or rate of change. In the context of a financial market, this external signal can be interpreted as 
an advertising, a fad or a public perception in favor or against one of the two possible trading 
behaviors, thus periodically breaking the symmetry of the system. 

We are interested in the modification of the social processes of opinion formation and prop- 
agation of information among the economic agents. So we are naturally led to introduce the 
external forcing in the social term of the transition rates, that is, in the herding coefficient h. 
Thus, we modify the transition rates ([T) as 

n:+(n,f) = 7z{n^n + l,t)^ {N - n) {a + h+{t)n) , 

(7) 

n {n,t) — %{n^n-l,t)= n{a + h^{t){N — n)) , 

where the herding coefficients are now different in the two possible directions and both are time 
dependent functions given by 

F 

h+{t)= ho + -cos{Q.t), 

(8) 

F 

h^{t)= /!o--cos(nO, 

with hf) playing the role of an average herding coefficient, F acting as the force or signal intensity 
applied to the whole system and Q. as the signal frequency. The reason for adding and subtracting 
the cosine term, respectively in h+{t) and h-{t), is just so that the herding coefficients are in 
anti-phase. Note that the forcing term is not proportional to the total force exerted on the system, 
but to the total force per agent, ^. The rationale behind this particular functional form of the 
signal intensity is basically a limited resources assumption: the resources spent in order to exert 
the force on the whole system are divided among its constituents, so that if the system size 
increases, the resources spent on each of the agents decrease. In order to keep the transition rates 
always positive, this intensity must satisfy ^ <ho. 

Proceeding in a similar manner as for the original model in section|2l and applying the same 
approximations, we find the new Fokker-Planck equation: 



dP{m,t) _ d 
dt dm 



( lam -F{\- OT^) cos(£2f ) ]P{m,t) 



1 

2 dm^ 



^ + 2ho{l-m^)jP{m,t) . 



Note that, compared to the previous Fokker-Planck equation (O, the herding coefficient h has 
been replaced by its average or constant part, Iiq. Furthermore, there is a new time dependent 
term inside the drift function, li{m,t) — —2am + F{l — m^) cos(i2f ). Again, the conventional 
transformation rule leads us, within the Ito form, to the Langevin equation describing the process. 



m = ^(m,f)+ ^jD{m)-^{t) 

[A'a ^^^^ 

= -2am + F(l -m2)cos(i2f)+ W +2/!o(l - m^) • ^ (f ), 

where ^ (f) is, as before, a Gaussian white noise. 

The constant part of the herding coefficient, Hq, plays exactly the same role as the herding 
coefficient itself in the original Kirman model. The new term of the drift function, on the contrary, 
changes fundamentally the general behavior of the system. In particular, it will try to force the 
system to follow the periodicity of the signal, as it will periodically favor or oppose the tendency 
towards m — Q caused by the first drift term. The new parameter F, the strength of the signal, 
modulates the intensity of this effect. Note also that the factor (1 — m^) inside the new drift term 
causes its effects to vanish at the extremes of the population configuration space and its absolute 
value to be maximal at its center for cos(i2f) = ±1. As a consequence, this new term does not 
help the system to exit the absorbing states, and thus some idiosyncratic behavior {a > 0) is still 
needed in order to observe transitions between both full agreement states. On the other hand, the 
equilibrium point of the drift term is no longer a constant at m — 0, but a periodic function of 
time around this central point. 

Concerning the competition between the deterministic and the stochastic terms of equation 
(|5]l, the inclusion of an additional external signal in (fTOl i has the general effect of periodically 
counteracting or enhancing the deterministic force. For a deeper understanding of the transition 
induced by the multiplicative noise upon the deterministic potential and the symmetry breaking 
role of the external signal, let us write an effective potential for the Fokker-Plank equation (|9]i, 

yeff(m,f) = (ho -a) ln(l -m^) -mFcos(i2f). (11) 



Again, we leave its derivation for [Appendix B Note that the new term, related with the input 



signal, is linear in the magnetization variable m, and thereby it breaks the symmetry (m o —m) 
of the potential, for cos(£2f) 0. 

The time evolution of the effective potential is shown in Figure |2] where five different time 
snapshots are presented for three values of the parameter a. For sufficiently small values of a 
the system is not be able to follow the signal, even if a double well effective potential has been 
induced, since these values are too small to allow for the system to exit the wells. Larger values 
of a but still smaller than ho allow both for the creation of a double well potential and transitions 
between them, while intermediate values of a larger than /iq give rise to a rather wide monostable 
effective potential where the system can still be largely driven by the external signal. Very large 
values of a lead to narrow effective potentials with minima moving closely around the center of 
the opinion space. 

The particular modification of the behavior of the system due to the introduction of the forc- 
ing signal depends on the specific values of the model parameters. Figure |3]contains four panels 
which show the effect of increasing the intensity of the signal from F = to its maximum al- 
lowed value F — NIiq. Note that the values chosen for the idiosyncratic switching and the herding 
parameters are, respectively, a = 0.001 and ho = 0.01, so we find a fairly bimodal phase, as in 
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Figure 2: Snapsliots of tlie effective potential V'eff(m,() in various time instants f = f* and for tliree different values of the 
idiosyncratic switching tendency, a = 10^^, 10^^, 10"' . The rest of the parameter values are Iiq = 0.01, F = 0.2 and 
Q, = 0.02. Note that the values of a shown here correspond to the three main cases a < ho, a = ho and a > hg. 

the unforced case. We observe the important differences between both upper panels (a and b) 
and those below them (c and d): the system does not follow the periodicity of the input signal 
in the former cases, while it certainly and quite accurately does in the latter ones. Obviously, no 
periodic behavior is expected in panel a, where there is no input signal (F ~ 0). Regarding panel 
b, characterized by a weak signal intensity, we observe a behavior quite similar to the previous 
unforced case, with no appreciable periodicity. Nonetheless, the particular trajectory of the sys- 
tem is clearly different and we notice a certain increase in the number of transitions between both 
extremal states, which point at the fact that a weak signal intensity is already able to facilitate 
some transitions. Even if in both panels c and d the system shows a clear order in time, it is 
also noticeable that in the first case, panel c, the periodicity is not so accurately followed as in 
the second one, panel d: some periods are longer or shorter than the others, some transitions are 
missing, there is a longer phase lag, etc. 

4. Resonance phenomena 

We present, in this section, the main results obtained from the analysis of the model intro- 
duced in section[3]conceming the existence of a resonance phenomenon, i.e., a maximum in the 
response of the system to the external signal as a function of a given control parameter As we 
have already pointed out, the capacity of the system to follow the input signal is basically related 
to the presence of some idiosyncratic behavior or random autonomous decision making in the 
system (a > 0). Therefore, we choose this idiosyncratic switching intensity a to be our control 
parameter, while we keep constant the average herding coefficient /iq and the forcing strength 
F. In particular, we fix their values as ho = 0.01 and F = 0.2, as in panel c of Figure |3] Fi- 
nally, we have chosen to show here the results for a number of agents — 200, although we 
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Figure 3: Black line: Time evolution of the opinion index simulated for different values of the signal intensity F and fixed 
values of the parameters a = 0.001, ho = 0.01, Q. = 0.02 and N = 200. Grey line: External periodic signal in arbitrary 
units (a.u.), intended only for period comparison. 



also performed simulations with N = 50 and = 800, finding a generally equivalent behavior as 
discussed below. The signal intensity per agent, ^, is thus only a 10% of the herding coefficient. 

Let us start by considering the influence that varying the idiosyncratic switching tendency a 
has over the trajectories or time series of the opinion index m (see Figure |4|i. Note that there is 
no maximum allowed value for this variable, as it was the case with the signal intensity F: its 
only constraint is that it must be a > 0. Therefore, we simply choose a reasonable range which 
includes the different behaviors described in the previous section and observed in Figure|2] from 
a fully bimodal case, with almost two deltas at the extremes of the probability distribution of 
states; up to a fully unimodal case, with an almost perfect Gaussian distribution. 

The first of these cases (panel a) corresponds to a value a = 0.0001 = j^. It is clear that, 
in this case, the system does not follow the periodicity of the input signal and, in fact, it stays 
pretty close to a full agreement state for most of the time. This is due to the extremely low 
level of idiosyncratic behavior compared to the herding intensity, which makes it very difficult 
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Figure 4: Black line: Time evolution of the opinion index simulated for different values of the idiosyncratic switching 
tendency a and fixed values of the parameters ho = 0.01, F = 0.2, Q. = 0.02 and N = 200. Grey line: External periodic 
signal in arbitrary units (a.u.), intended only for period comparison. 



for the system to exit these consensual states and, therefore, the transitions between both of 
them cannot be as frequent as the forcing tries to impose. Even if the example shown in panel 
b, corresponding to a value a = 0.001 = y^, was as well in the bimodal region in the original 
model, the idiosyncratic behavior intensity is already strong enough so as to take the system out 
of the absorbing states with a frequency allowing it to fit the external signal in almost all the 
transitions. Note, for instance, that the last switch is missing in the trajectory. This is not the 
case in panel c, where the idiosyncratic coefficient equals the herding intensity, a ~ 0.01 — ho, 
and thus the probability distribution of states is uniform in the unforced model. We observe 
here that the system easily follows the periodicity of the input signal, although the trajectory 
starts to be fairly noisy. The case shown in panel d, with a = 0.5 — 50 • ho, corresponds to a 
rather unimodal probability distribution of states in the original model. The output trajectory 
accurately follows the period of the forcing signal. However, the great strength of the first term 
of the drift function, the one proportional to a (see equation (fTOli). prevents the system from 
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collapsing to the consensual states. Therefore, the competition between the first and the second 
terms of the drift function keeps the system oscillating in the proximity of the balanced situation, 
m — 0. 

We note that, by applying a driving signal whose absolute strength per agent is only a 10% 
of the herding intensity, we can obtain (for instance in panel c of Figure |4|i a large amplification 
of this input, as the system tends to collapse to the extremal states m = ±1. This is because, 
on the one hand, the driving signal breaks locally in time the symmetry of the problem, by 
favoring periodically one or the other states, and on the other hand, the stochastic or diffusive 
term (see equation ( fTOl l) tends to make any asymmetry grow towards a total consensus, as has 
been reported above. However, this effect is not observed for every value of the idiosyncratic 
switching parameter a: it is clearly not present in panels a, where there is even no periodic 
behavior, and d, where there is no collapse to the extremes. In short, what we observe is a 
maximum in the amplification of an external driving signal as a function of the intensity of the 
agents' idiosyncratic behavior, a phenomenon which correctly fits the definition of resonance. 

In order to quantify this resonance phenomenon, we consider the quality of the global re- 
sponse following the input forcing, that is: how much the system output m fits the input signal 



jj cos (Q.t). For this purpose, we use the spectral power amplification factor r\ dGammaitoni et al 



mi, 



^ = -\{e-''^m{t))\\ (12) 

which is basically the normalized amplitude of the Fourier component at the forcing frequency 
£2 of the time series m{t). Note that the signal intensity used for the normalization is the absolute 
value of the total force applied to the whole system, F . A larger T] is related with a better 
entrainment of the system by the external signal. Note that a larger entrainment can be caused 
both by a better fit of the input signal or by its amplification. 

The results obtained for the spectral power amplification factor Tj are shown in Figure |5] for 
three different system sizes. The same general behavior is observed in all of them: there is a 
region in the parameter space in which the system response to the weak periodic driving signal 
displays a maximum value as a function of the idiosyncratic switching intensity a. Furthermore, 
the particular values of these maxima are quite large, all of them implying a large amplification 
of the input signal. Note as well that this amplification occurs for a fairly wide range of values 
of the control parameter a in all three cases. These results suggest that there is a certain range 
of idiosyncratic behavior intensity which largely improves the entrainment of the system by an 
external periodic signal. 

The only difference worthy of remark between the three examples in Figure|5]is due to the 
granularity term, the only one dependent on A^. In fact, the three sizes collapse in the same curve 
for higher values of a, while they are clearly different for intermediate values and join again 
together for the lower a region. This behavior is caused by the unequal importance of the two 
a dependent terms in the Langevin equation (fTOl i: the granularity and the first drift term. For 
very small values of the control parameter a, the system stays most of the time absorbed in the 
consensual states, so neither of these terms has any noticeable importance. In the intermediate a 
region, the granular dependent term plays a relevant role in taking the system out of the con- 
sensual states with a probability uncorrected with the frequency of the input signal. Therefore, a 
large granularity term — small number of agents — leads to a worse amplification of the input 
signal and, thereby, to a lower spectral power amplification factor r]. Finally, for large values of 
the control parameter, the system does not even reach the absorbing states and the behavior is 
predominantly led by the first term of the drift function, so the granular effects become negligi- 
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ble. Note as well that the relative difference between the curves is much smaller when comparing 
the = 800 with the = 200 cases than when comparing this latter with the = 50 example. 
This is simply a consequence of the granularity term being inversely proportional to the system 
size A^, so for larger and larger A^ the results become more and more similar. 

Let us now discuss some insights regarding the particular features of the resonance mecha- 
nism f ound. Indeed, this mechanism is reminiscent of a ph enomenon known as stochastic reso- 



nance (iBenzi et all 1 1 98 IL 1 1 982t iNicolis and NicolisL 1 1 98 11) . since the maximum in the response 



of the system to the external signal is related to the relative importance of the stochastic term 
as compared to the deterministic one. Note, however, that in our case (see equation (fTOl l) this is 
only true for A^ tending to infinity. Furthermore, if in a stochastic resonance the control parameter 
would normally be the intensity of the stochasticity (related to h in our case), we have chosen 
the intensity of the deterministic term, a, as our control parameter Both the stochastic resonance 
and the one here presented take place when a deterministic potential is driven by a combination 
of a periodic and a random forcing. Nevertheless, the roles played by this random force and the 
deterministic potential are fundamentally different in our case: both factors exchange their roles. 



In the case most commonly found in the literature (iGammaitoni et al.Lll998i) . the potential has a 



double-well shape and the noise intensity or diffusion function is additive and is a constant for 
the whole variable range. On the contrary, the case here presented is characterized by a monos- 
table deterministic potential and a multiplicative noise intensity which is maximum at the center 
of the opinion space and vanishes at its extremes. Therefore, the deterministic force pulls the 
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system towards the center while the net effect of the random forces is to drive the system out of 
the center and towards the extremes, in contrast to the usual stochastic resonance scenario. 



5. Conclusions 

The aim of the present paper is to advance towards a quantitative understanding of the in- 
fluence of an external source of information upon a financial market characterized by a certain 
herding behavior The stochastic formalism used as a point of departure for our investigation 
incorporates individual behavioral heterogeneity as well as a tendency for social interaction, in 
the tradition of Kirman's seminal ant colony model. We have extended this original model by 
adding a periodical signal modeling a varying incoming information from an external source. 

A transition takes place in the original Kirman model from a monostable to a bistable be- 
havior when increasing the herding propensity of the agents with respect to their idiosyncratic 
tendency. We have reinterpreted this noise-induced transition in terms of the mono- or bistability 
of an effective potential. The monostable case can be understood as a market where each of both 
possible strategies is used by approximately half of the traders, while the bistable configuration 
corresponds to a market where there is always a clear majority of traders using a given strategy, 
even if the chosen one can change over time. In this context, we have shown that, when an input 
of external information is taken into account, a resonance phenomenon may take place under 
certain circumstances. In particular, we have found a certain range of values of the idiosyncratic 
behavior intensity a that maximizes the response of the system to a weak periodic input signal. 
The specific range of values for which an amplification is observed depends on the intensity and 
the frequency or rate of change of this incoming information. 



Within a financial market framework like the one presented bv lAIfarano et al.l(l2008h . a max- 
imal amplification of the external signal would simply mean that the incoming information influ- 
ences virtually the whole market and is able to convince most of the traders to a particular mood, 
to act in a certain manner or to use a specific trading strategy, whether it is an optimistic or pes- 
simistic mood, a buying or selling activity, a chartist or fundamentalist strategy, etc. We have 
shown the existence of three different regimes regarding the capability of the system to follow 
the external signal. In financial terms, we could translate these regimes as: 

1 . Impossibility for the market to assimilate the incoming information in a reasonable time, 
that is, in a time scale shorter than that of the variations of the information itself (values of 
the idiosyncratic switching intensity a below the range of the resonance). 

2. Exaggeration of incoming shocks of information, i.e., even small outside changes lead to 
large market movements (values of a within the range of the resonance). 

3. Precise assimilation of incoming information but with a very small effect in the market 
(values of a above the range of the resonance). 

The amplification of small external signals can lead to enormous price variations, a large volatil- 
ity and, therefore, to a great instability of the market. The results presented in this paper support 
the idea that the greater intensity of the herding with respect to the idiosyncratic tendency may 
play an important role in the development of such instabilities in a financial market open to 
external incoming information. 
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Appendix A. Simulation method 

In contrast with the methods present ed in the previou s literature, we have used a Gillespie 



algorithm for the simulation of the model (lGillespieLll977h . A single realization of this algorithm 



represents a random walk with the exact probability distribution of the master equation, therefore 
generating statistically correct trajectories of the stochastic equation. In particular, we find the 
unbiased times when the transitions take place. 

For a stochastic system with time dependent transition rate^ ^^(0 = « ± 1,0' '^^^ 

given that the last transition took place at time fi, the cumulative probability of observing an 
event « — n ± 1 at time t2 can be written as 

and we can equate this expression to a uniformly distributed random number between and 
1 in order to find an equation for the time of the next event, tj, 

'2 



n^{t)dt = -ln(l -M*) = -ln(M*), 
t\ 

where we have used that and 1 — are statistically equivalent. By that means, two different 
times are found: , corresponding to the transition rate K'^{t), and fj^, related to the transition 
rate ?r^(f)- The event actually taking place is the one with a shortest time. 

Applying the last result to the transition rates of the original Kirman model, as given in 
equation ([T]i, and solving for t2 we find 

, ln(M+) 



[N -n) (a + hn)' 



n {a + h{N -n)) 

Regarding the forced case, whose transition rates are given by (|7]) and we find 

, InfM^) F n / ^ , Ls\ 

h = fi - 7 r-r-^ r + (sm(nfi) - sm(£2f7)) , 

2 {N~n){a + hon) NO. a + hon^ ^ ' 2 1)^ 

\\\{u-) F (N-n) , , , , 
U = fr ^—-^ — H — - — - — - (sm(i2f, ) - sm(i2fi )) , 

2 ' n{a + hoiN-n)) ND. a + ho{N-n)^ ^ ^' ^ i;;, 

equations whose solution can be approximated by different numerical methods. 



''These rates may also depend on other variables, as it is the case in the Kirman model, but we are here interested only 
in the time dependence. 
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Appendix B. Effective potential 



We derive, in this appendix, the effective potential both for the original and the forced Kir- 
man dynamics, presented in the equations (|6]l and (fTTT l. Let us start by defining the "effective 
potential" V'eff('«) as 



Pst{m) = ^"'exp 



-1 / Veff(w) 



D 



where Pst('w) is the stationary state probability distribution, D is an effective noise intensity and 
the constant plays the role of a normalization factor. 
For the general Fokker-Planck equation 

1 ' = - ^ [q{m)P{m,t)] + ^ [Dg{m)P{m,t)] , 
at dm oni'^ 

the stationary distribution is found by assuming -gf- — and solving the resultmg equation. By 
this means, the general effective potential can be written as 

. - / + 2 I '-f^^,dm, (B.l) 

J g{my J dm g[m) 

and, applying this definition to the Fokker-Planck equation (|3j, the particular effective potential 
for the Kirman dynamics is found to be 

Veff(m) = {h - a) ln(l - m^). 

Note that this effective potential Veff (m) is not to be confused with the deterministic potential, 
which is always monostable and can be found by simply integrating with respect to m the deter- 
ministic part of equation (|5]l. 

Even though in the case with an external periodic forcing it is not possible to write a stationary 
state probability distribution, we assume that, at any point in time, the decay of the system to a 
"quasi-stationary state" is faster than the variation of the input signal. Therefore, we keep the 
previous definition of the effective potential as an approximation to this time dependent case. 



P{m,t) « ^"'exp 



D 



Thereby, applying equation (IB. lb to the Fokker-Planck equation (|9|l leads to the particular func- 
tional form 

Veff(m,f ) = {ho — a) ln(l — m^) — mF cos{Q.t) 
for the forced Kirman dynamics. 
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